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This  report  .■as  prepared  by  the  Department  of  Civil  Engineering  , 
Cniv-.rsity  of  Florida,  Gainesville  FL  32511,  under  Contract  Number 
'“'08h35 -83-C-0136 ,  Task  No.  83-1,  for  the  Air  Force  Engineering  and  Services 
C  ..  .er,  Engineering  and  Services  Laboratory  (HQ  AFESC/RD) ,  Tyndall  Air  Force 
.  ,  Florida.  This  report  is  published  as  submitted  to  the  University  of 
I  iorida  by  Captain  Dale  M.  Bradley,  as  his  Master  of  Engineering  Tb.esis, 

...nd  r  the  direct  supervision  of  Professors  F.  C.  Townsend,  F.  E.  Fagundo, 
and  J.  L.  Davidson.  Captain  Paul  L.  Rosengren,  Jr.,  was  the  HQ  AFESC/RDCS 
‘.tOje.'.r.  Officer. 

This  report  surnisrines  work  performed  between  July  and  December  1983  and 
discus sen  the  feasibility  of,  and  scaling  relationships  for,  centrifugal 
.cocoil  br.g  of  the  GLCM  shelter.  Similitude  between  model  and  prototype  is 
deemed  possible  if  constitutively  similar  reinforced  concrete  and  soil  are 
used  in  the  model.  However,  currently  available  US  centrifuges  limit  GLCM 
shelter  model  scales  to  1/150  to  1/3J0,  which  caused  difficulties  in  model 
construction  with  microconcrete  and  miniaturized  reinforcement.  Therefore, 
centrifugal  r.odeling  of  GLCM  components,  coupled  with  computer  analyses, 
is  recommended. 

This  report  has  been  reviewed  by  the  Public  Affairs  Office  (PA)  and  is 
rr-leasable  to  the  National  Technical  Information  Service  (NTTS)  .  At  NTIS 
it  will  be  available  to  the  general  public,  including  foreign  nationals. 


This  technical  report  has  been  reviewed  and  is  approved  for  publication. 
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SECTION  I 


INTRODUCTION 


A.  GENERAL 

Appropriate  load  and  response  testing  is  fundamental  to  design  of  a 
new  structure.  This  type  of  evaluation  is  done  commercially  with 
component  testing  and  analytical  computer  simulation.  A  much  more 
intensive  and  comprehensive  testing  program  is  needed  for  protective 
structures  developed  for  the  Department  of  Defense.  The  degree  of 
uncertainty  must  be  quite  narrow  to  satisfy  defense  design  engineers. 
Simple  application  of  safety  factors  is  not  sufficient  and  the  entire 
structure  should  be  tested  under  explosive  loading.  This  often  requires 
full-scale  testing  of  structures  over  several  geological  sites  to 
accurately  evaluate  the  design  assumptions  and  computer  models.  This 
has  been  done  in  the  past,  but  at  monumental  cost  and  substantial  safety 
risk. 

To  reduce  testing  costs  and  safety  risks,  the  Department  of  Defense 
decided  to  test  scaled  models  on  the  order  of  one-half  to  one-twentieth 
the  size  of  the  prototypes.  While  this  order  of  scaled  testing  reduced 
costs,  the  reduction  has  not  been  enough  to  meet  budget  criteria. 
Additionally,  there  are  some  technical  concerns  regarding  testing 
smaller  than  one-tenth  the  prototype  size.  Murphy  (Reference  1}  suggests 
that  the  dead  load  effects  of  the  structure  itself  are  not  appropriately 
handled  in  one  gravity  testing.  Since  dead  loads  play  a  substantial 
role  in  dynamic  response,  this  type  of  scaled  testing  is  limited. 
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Facea  wilh  these  concerns,  the  Air  Force  has  started  to  investigate 
the  possibility  of  using  a  centrifuge  to  test  its  new  Ground  Launch 
Cruise  Missile  (GLCM)  shelter,  (Figure  1).  This  structure  is  278  feet 
long,  164  feet  wide  and  38  feet  tall.  With  present  centrifuges  the  GLCM 
could  be  tested  at  1/150  to  1/300  the  size  of  the  prototype  with 
explosives  reduced  by  1/n  ,  where  n  is  the  number  of  gravities 
(Reference  2).  Coupled  with  data  provided  by  Mikasa  and  Takada 
(Reference  3)  which  indicate  that  the  centr’fuge  appropriately  models 
dead  load  effects,  the  centrifuge  shows  substantial  potential  for  cost 
reduction.  However,  major  concerns  prevail  in  the  community  regarding 
the  use  of  the  centrifuge. 

These  major  concerns  are  constructability  of  such  a  small  model  and 
applicability  of  model  data  to  prototype  response.  This  report  will 
evaluate  the  scaling  laws  that  will  ensure  the  similitude  between  monel 
and  prototype.  To  the  extent  possible,  the  formulated  scaling  laws  will 
indicate  the  feasibility  of  the  centrifuge  for  GLCM  shelter  testing. 

■3.  OBJECTIVE 

The  primary  objective  of  this  effort  is  to  investigate  the 
feasibility  of  using  a  centrifuge  to  test  scaled  models  of  the  GLCM 
shelter.  The  basic  premise  of  the  study  is  to  see  if  a  scaled  GLCM 
shelter  can  be  tested  in  a  centrifuge  in  such  a  manner  that  ensures 
similitude  between  model  and  prototype.  In  other  words,  can  the 
generated  data  of  a  centrifuge  model  testing  program  be  used  to  predict 
the  GLCM  shelter  response.  To  meet  this  objective,  the  following  must 
be  completed: 

1.  All  relevant  parameters  of  the  structure  and  surrounding 

material  that  affect  structure  response  due  to  explosive  loading 
must  be  identified. 
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2.  Using  the  response  parameters,  formulate  scaling  laws  that  will 
re1 ate  the  model  to  the  prototype. 

3.  Based  on  the  scaling  laws,  determine  the  physical  limitations  of 
the  model  design  and  testing  procedure. 

4.  Based  on  physical  limitations,  recommend  the  model  scale  and 
centrifuge  to  be  used  in  a  testing  program. 

K.  SCCPE 

Tne  effort  began  with  a  review  and  study  of  theory  of  model  ling 
techniques  and  procedures.  The  review  included  simple  to  complex 
problems  '■esult'ng  in  the  identification  and  use  of  the  Buckingham  Pi 
Theory.  The  theory's  application  is  .discussed  and  demonstrated  in  a 
generic  pronlen  and  used  to  prove  a  common  table  of  scaling  laws.  This 
established  the  methodology  to  develop  tne  uLCM  centrifugal  scaling 
laws.  The  response  parameters  were  then  identified  and  used  to 
establish  the  system  (independent)  parameters  wnich  would  affect  them 
due  to  the  leading  type  and  materials  involved.  The  research  progressed 
into  a  trial  and  e-ror  procedure,  using  the  Buckingham  Pi  Theory  to 
identify  a  complete  set  of  scaling  laws  to  govern  the  model  testings. 
Based  on  the  literature,  review  qualifications  and  limitations  to  these 
laws  are  addressed  and  clarified  based  on  the  extent  that  they  may  cause 
perturbations  to  similitude.  A  comparison  of  tne  available  centrifuges 
in  the  United  States  to  the  ohysica!  size  of  the  GLCM  shelter  was  then 
conducted  to  evaluate  possible  scales  to  be  used.  This  resulted  in  a 
list  uf  ootenf.ial  scales  and  centrifuges  to  be  matched  together  in  a 
model  testing  program.  The  effort  concludes  with  suggestions  for 
explosives,  instrumentation,  and  materials  to  be  used  in  the  model. 
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D.  BACKGROUND 


Practically  all  the  literature  dealing  with  scaling  uses  the 
Buckingham  Pi  Theory  for  the  basis  of  their  study,  regardless  if  it  is 
in  a  one-  or  multiple-gravity  environment.  A  substantial  list  of 
references  deals  with  scaling  laws  associated  with  one-gravity  tests 
(see  Reference  List);  however,  dynamics  centrifuge  research  is 
practically  nonexistent.  Several  researchers  (References  3, A, 5)  have 
studied  structures  over  soil  in  multiple-gravity  environments  but  only 
under  static  loads.  Many  others  have  studied  dynamic  loading  on 
structures  over  soil  but  in  one  gravity  (References  6,7).  Apparently 
only  a  very  few  have  studied  the  scaling  laws  in  the  condition  of 
dynamic  structure/soil  interaction  within  a  centrifuge  (References  8,9) 
Regardless  of  the  type  of  experiments  researched  or  conducted,  tne 
primary  axiom  of  scaled  testing  is  to  use,  when  possible,  the  same 
material  properties  in  the  model  as  in  the  prototype.  This  drastically 
simplifies  the  development  of  the  scaling  laws.  The  axiom  has  resulted 
in  the  development  of  both  microconcrete  and  gypsum  mortar  to  satisfy 
scaled-down  concrete  mixes.  For  soils,  simple  reduction  in  particle 
size  has  proven  unsatisfactory,  due  to  the  complex  interrelationship  of 
soil  strength  and  response  to  particle  size.  Subsequently,  researchers 
have  used  the  exact  same  soil  in  model  as  in  prototype  to  ensure 
consistency  in  material  properties.  Perturbations  from  this  action  hav? 
proven  insignificant  if  the  grain  size  is  maintained  much  smaller  than 
the  scaled  structure  element  (.Reference  8).  Previous  researchers 
(Reference  9)  feel  that  creditable  scaling  laws  can  be  developed.  This 
effort  will  concentrate  on  that  endeavor,  but  the  author  premises  such 
research  with  the  questionable  constructability  of  the  complex  GICM 
shelter  at  any  scale  smaller  than  1/100. 


5 


SECTION  II 
MODELLING  THEORY 


A.  GENERAL 

In  the  initial  evaluation  ot  models,  it  was  necessary  to  define  "he 
function  that  relates  all  the  parameters  affecting  the  predicted 
phenomena  (dependent  parameter).  If  stress  (q j )  is  dependent  on  q2,  q2, 
qq  (system  parameters)  then  there  exists  a  function  such  that: 

Oj  =  03*  94)  (]  ) 

In  simple  systems  in  which  the  dependent  parameter  (qj)  is  only 
dependent  on  a  small  group  of  parameters  (q2,  q2,  qq),  the  function  is 
readiiy  evaluated;  however,  in  complex  systems,  such  as  the  GLCM,  the 
list  of  parameters  is  quite  numerous,  making  the  determination  of  the 
function  practically  impossible. 

In  1914,  E.  Buckingham  proposed  a  theory  (known  today  as  Buckingham 
Pi  Theory)  which  eliminates  the  neea  to  know  the  function  and  allows 
direct  similitude  through  the  use  of  *  terms.  Buckingham  states  that  if 
the  equation  relating  the  dependent  parameter  to  system  parameters  is  a 
dimensional ly  homogeneous*  it  can  be  reduced  to  a  complete  set  of 

dimensionless  products  (Reference  10).  Therefore,  if 

Ql  s  f(q2.  Q3»  Q4»  •  •  •  qn)  (2) 


r, 

2 


a 

v.  $ 


*  To  be  dimensionally  homogeneous  is  to  be  completely  defined  by  one 
set  of  relevant  dimensions  such  as  mass,  length,  and  time. 
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it  can  re  defined  as 

'(d'l  c .}  i  q  - ,  q  i ,  •  •  *  y  y  ~  0 

For  wh  ch  the  Buckingham  Pi  Theo'y  rtates, 

"If  the  equation  f(q, ,  q-,  q3,  q4...qnJ  =  0  H  complete. 


the  solution  nas  *  rorm  *9,  tt?. 


.  tt  .  )  =  0  wh<ve 

n  -  ;< 


the  it  terms  a  .-  i raereroenr  produ “t S  or  toe  iMremeters 
q,,  q?,  etc.,  l  id  are  <1 1  mens  i  r.n !  es?  ir,  t-w  ,rur  :o  recta 1 
dimensions"  (P  ;fermc?  '.!).  In  this  case.  '<  t  the  r-nk  of 
tne  diner si  or/ 1  matrix. 

tssest'.ii  iy ,  Hue*  i  rgiam  is  saying  that  if  all  the  re  i  event 
parameters  are  included  in  tne  initial  equation,  then  si  nil  it  tide  tot.-  ".-.i 
model  ana  protetype  can  be  achieved  if 

-r  =  t  Mi 

n  p  ' 

"model  p^ototyoe 

b.  FX AMPLE 

Langhaar  (Reference  10'  provides  an  illustrative  example  of  tne  use 
of  the  Bucxingham  Pi  Theorem  on  a  hypothetical  problem  in  which  too 
dependent  parameter  to  be  evaluated  is  P.  P  is  dependent  on  Q,  R.  S.  T, 
U,  &  V  (System  Parameters}.  Wnere  p,  0,  P,  S,  T,  U,  5  V  are  defined  in 
terms  of  mass  (M),  length  (L),  and  time  ft)  as 


p  =  (MZl)\  Q  -  R  =  (£)\  S  =  (t3)\ 


jKr  ,K-  />  o  K-j 

I  =  (lA)  5,  U  =  (4-1  and  V  =  (ML“tc) 


M“t 


I 


/ 


lEt'Waf 


m  jnuwni 


f  i 
* 

/  I 


With  all  parameters  defined  in  fundamental  units  M,  l,  t,  a  dimer si  or: i 


matrix  is  developed. 


M 


To  e n-urc*  tnat  the  '>  terms  are  1  ndependenf. .  the  rcws  of  the  dimensional 
matrix  must  be  linearly  independent.  If  the  rank  (r)  of  the  matrix  is 
usual  to  the  number  of  rows,  independence  is  guaranteed.  Tne  size  of 
any  internal  square  matrix,  whose  determinant  is  nonzero,  indicates  mat 
rank.  This  value  determines  the  number  of  »  terms  (s)  to  he  generated 
wher° 


?  R  S  T 


2 

-1 

3 

0 

0 

-2 

1 

* 

1 

c 

-i 

0 

2 

1 

2 

0 

i 

0 

0 

1 

_  1 

2 

and  a  equals  one  number  of  parameters  in  the  problem.  The  determinant 
of  the  last  three  co'um.ns  (T,  U,  V)  is 


0 


-2  1 

2  1  ? 

1  -I  2 


=  1 


Therefore , 

1}  The  rank  equals  3. 

2)  The  rows  are  linearly  independent. 

3)  The  number  of  expected  *  terms  is  4. 


,  p  ) 


Since  the  «  terms  are  products  of  the  parameters  they  have  the  form  of 

K,  K,  K,  K.  Kc  K7 
it  =  P  1  Q  2  R  3  S  4  T  5  U  6  V  7  (7) 
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and  since  they  must  be  dimensionless,  the  following  must  be  true: 


M°L°t°  =  (M?t)Kl^)K2(^)K3(t3)K4(L2t)Ks(4-)K6(ML2t2)K7  (8) 

L  Me 


By  combining  units  the  following  is  true: 


0,  0  0  (2KrK2+3K3-2K6+  Ky)  (K,  “K3+-K5  +K6  +2K/)  (K2  +3K4+K5  ~Ks  +2M  ^ 


2K1-K?+3K3-2K6+K? 

=  0 

(10) 

K1-K3+2K5+K6+2K7 

=  0 

(11) 

K2+3K4+K5-K6+2K7 

=  0 

(1?) 

which  are  the  homogeneous  linear  algebraic  equations  whose  coefficients 
are  the  numbers  in  the  rows  of  the  dimensional  matrix. 

Now  choose  the  repeating  variables  equal  in  number  to  the  number  of 
rows  in  the  dimensional  matrix.  These  variables  are  system  parameters 
not  including  the  dependent  parameters.  They  should  meet  the  following: 

1.  independent  within  themselves 

2.  easily  controlled  experimental ’y 

3.  contain  all  fundamental  units  at  least  once 

4.  not  equal 


langhaar  (Reference  10)  chose  Kg, Kg,  and  Kj.  Solve  the  algebraic 
equations  for  the  exponents  of  the  chosen  repeating  variables  —  Kg,  Kg, 
and  K 7.  Starting  with 
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resulting  in 


K5  =  -llKj  +  9K2  -  9K3  +  15K4  (15) 
K6  =  5Kj  -  4K2  +  5K3  -  6K4  (16) 
K?  *  8Kj  -  7K2  +  7K3  -12K4  (17) 

Then  assume  values  for  ,  K2,  K3,  and  K4  such  that  only  one 
appears  in  each  it  term.  This  results  in  the  Matrix  of  Solutions. 
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This  states  that  the  prototype  response  Pp  relates  to  the  model  response 
pm  as 


P 

P 


U  5V  8T  11 
(  m  m  P  Id 
Pm 

P  'P  in 


( 26) 


In  the  same  manner,  all  other  model -to -prototype  relationships  can  be 
determined. 


SECTION  III 


STANDARD  SCALING  DERIVATION 


A.  DISCUSSION 

The  vast  majority  of  the  work  cited  in  the  literature  assumes  that 
the  sane  materials  will  be  used  in  the  model  as  in  the  prototype  or  that 
at  least  the  material  properties  will  be  kept  constant.  Contingent  on 
this  assumption,  a  standard  scaling  table  is  o';ten  presented  but  never 
derived.  A  sound  understanding  of  these  standard  scaling  laws  is 
required  for  any  modelling  study.  It  is  therefore  beneficial  to  derive 
and  evaluate  these  laws  and  not  take  them  at  face  value  as  done  in  other 
model  studies.  A  common  form  of  this  scaling  table  is  Table  1.  It  is 
essentially  the  same  table  presented  in  Nielson's  (Reference  9)  and 
Cheney's  (Reference  12)  work.  It  lists  the  relevant  physical  quantities 
from  length  to  hydrodynamic  time.  All  must  be  appropri ately  defined  in 
a  model  study. 

Material  strength  parameters  are  omitted  because  they  must  be 
constant  between  model  and  prototype  for  the  table  to  be  valid.  Each  of 
the  scaling  laws  corresponding  to  its  respective  physical  quantity  is 
derived  from  an  appropriate  dimensionless  product  (n  term).  If 
all  »  terms  were  presented  for  each  of  these  quantities,  it  would  be 
apparent  that  material  properties  strongly  influence  the  terms. 
Subsequently,  researchers  such  as  Sabnis  (Reference  13),  Morris 
(Reference  8),  Nielson  (Reference  9),  and  many  others  fix  these 
properties  and  keep  them  constant  in  model  and  prototype.  This  causes 
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TABLE  1.  STANDARD  SCALING  LAWS 


QUANTITY 

LINEAR  DIMENSION,  (L) 

GRAVITY  (g) 

AREA,  (A) 

VOLUME,  (V) 

DYNAMIC  TIME,  (trf) 

VELOCITY  (DISTANCE/TIME),  (v) 
ACCELERATION  (DISTANCE/TIME2),  (a) 
DENSITY  (HASS/VOLUME ) .  (p) 

UNIT  WEIGHT  (FORCE/UNIT  VOLUME)  (r) 
FORCE,  (F) 

STRESS  (FORCE/AREA),  (a) 

MASS,  (M) 

ENERGY,  (E„) 

STRAIN  (DISPLACEMENT/UNIT  LENGTH),  (c) 
HYDRODYNAMIC  TIME,  (th) 

IMPULSE 


FULL  SCALE 
(PROTOTYPE) 

I 

1 

1 

1 

I 

I 

1 

I 

1 

1 

1 

1 

1 

1 

1 

1 


CENTRIFUGE 
MODEL  AT  n  q‘s 

i/n 

n 

1/n2 
1/n3 
I  In 

1 

n 

1 

n 

1/n2 

1 

1/n3 

1/n3 

1 

1/n2 

1/n3 
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the  scale  of  length  and  scale  of  gravity  to  be  the  same.  A;  Section 
III  h  indicates,  this  allows  an  accurace  and  simple  formulation  of 
scaling  laws  listed.  Note  that  this  derivation  does  not  use  it  terms 
exclusively.  Only  the  length  scale  and  the  gravity  it  term  are  needed. 
The  remaining  quantities  fall  out  from  common  definition. 

B.  TABLE  i  DERIVATION 

This  section  lists  a  step-by-step  formulation  of  Table  1.  Each 
scaling  law  derived  is  formulated  based  on  previous  laws  or  defini¬ 
tions.  Notp  that  the  subscripts  m  and  p  correspond  to  model  and 
prototype,  respectively. 


1.  LENGTH  (L) 

By  definition  the  scale  factor,  =  j-E- 

m 


(27) 


2.  GRAVITY  (g) 

By  using  a  *  term  which  includes  L  and  g,  such  as  - 4i---  where 
E  is  modulus  of  elasticity,  the  following  is  true: 


Lap 
nvm  m 


L  q  p 
PTP 
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By  assuming  the  same  materials  —  E  =  F 

m 


p  and  pm  =  pp  the  following 


is  true: 


9m  =  •  r- =  s 


\j  L  L 

m  m 


9m  =  SLgp  (28) 


•  •  =  n,  where  n  *  #  of  gravities. 

3.  AREA  (A) 


(29) 


A  =  L 
m  m 


L  =  J. 
m  n 


L  A 

a  =  4-  =  4 

m  T 

n  n 


\  =  "jf  (30) 


4.  VOLUME  (V) 


V„  '  L/  •  L„  ■/.»!>  •  LP3 


Ln  V 

v  =  P  =  P 

m  ^  4 


5.  DYNAMIC  TIME  (td) 


“m'-j  (3D 

n 


L  -  1/2  g  t 
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9 


Lm  3  ?r  •  \  -  n9p 
2L 


P  t 


-JL.  =  _£ 

nyp  n 


t  -  — £ 
dm  n 


(32) 
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/  . 


i 


6.  VELOCITY  (v) 

L 


v  '  um/t  ,  L  =  _£.  t  =  JL 

rn  m  n  *  ~  — 


m 


n  m  n  m  n 


Lp/n/t  .  -  -£  =  v 

p/n  tp  vp 


vm  =  vp  (33) 


7.  ACCELERATION  (a) 


m  ^m/t^2  •  in  n 


L  t 

I  =  P  t  ~  P 

*  in  7T~  »  ~  ~ 


m  n 


a  =  p/n/t  2,  2  «  iig.  •- 


n 


p/n  F7 


na 


am  =  n3n  (34) 


8.  DENSITY  (p) 


MASS/UNIT  VOLUME 

This  is  constant  in  all  gravity  environments. 


P  s  D 

m 


(35) 


9.  UNIT  WEIGHT  (y)  FORCE/UNIT 


VOLUME 


Ym  =  Vm  .  g„  *  ngp  ,  pm  =  p 


Ym  -  ng  p  =  ny 

m  3p  p  '  p 


Ym  ”  %  (36) 
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10.  FORCE  (F) 


F 


m 


Y  V 
m  m 


F 

01 


nr  V  ,  . 
P  P/n3 


(37) 


11.  STRESS  {a) 


(38) 


12.  MASS  (M),  with  respect  to  dead  load 


M 

m 


F 


F 

m 


(39) 


13.  ENERGY  (tn) 


19 


16.  HYDRODYNAMIC  TIME 

T  (H  2) 

+  mv  m  ' 

V  ^ - 

m  V 

m 

(th) 

Assuming  sane  material: 

a>  Tm  3  Tp 

Time  Factor 

b)  Cv  =  Cy 

m  yp 

Coefficient  of  Consolidation 

c)  Hm  =  »p/ 
y‘  n 

Drainage  Path 

(43) 


•  t  . 


T 

_  P 


H  2 
P 

7 


t 

h 


tH  -  -y-  (44) 

m  n 


17.  IMPULSE  z 

' r'"  tan  ■  ^  •  $  ■  \  ■  -r 


F  _  d  I 
I  =  P  _P  ,  _E 
m  72  n  3 
n  n 


(45) 
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SECTION  IV 
PISETS  PROGRAM 

To  automate  the  *  term  generation  procedure  outlined  in  Section  II, 
Theodore  Self,  a  graduate  student  at  the  University  of  Florida,  wrote  c 
API  computer  program  called  PISETS.  This  program  was  written  to  run  or 
the  Commodore  Super  PET  mi crocomputer,  usi ng  Waterloo  Computing  Systems 
micro  API  version  1,  and  on  the  Northeast  Regional  Data  Center  system 
through  the  Virtual  Machine/System  Product  Conversion  Monitor  System 
using  VS  API.  A  computer  user's  manual  and  program  listing  are 
contained  in  Appendix. 

The  program  requires  the  dimensional  matrix  and  fundamental  units  < 
input.  Using  the  procedure  presented  in  Section  II  it  generates  the 
solution  matrix  for  all  or  selected  tt  sets.  With  system  parameters 
numbering  10  or  more  the  generation  of  the  solution  matrix  takes  a 
considerable  effort.  Noting  that  each  set  of  independent  repeating 
variables  generates  one  *  set  and  realizing  that  there  are  the  number  c 
system  parameters  factorial  possible  *  sets,  the  use  of  PISETS  is 
imperative. 

Use  of  the  program  is  limited  only  to  the  available  computer 
internal  working  memory  (RAM).  Processing  time  is  just  a  few  seconds. 
Printed  output  is  minimal  unless  a  large  dimensional  matrix  is  used  am 
all  possible  t  sets  are  needed.  Unfortuately ,  the  standard  SuperPET 
printer  is  not  compatible  with  APL.  The  mainframe  must  be  used  to 
obtain  a  hard  copy.  Presently  a  3  by  22  matrix  has  been  run  on  the 
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above  systems  with  no  difficulty  Nnfp  that-  ™ 

ly*  Note  that  more  than  the  standard 

three  fundamental  units  can  >e  used  if  th D  i 

used.  IF  thermal  or  electrical  units  are 

needed  in  the  study.  PI  SETS  will  accomodate. 
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SECTION  V 


INITIAL  PROBLEM  SET-UP 

A.  IDENTIFICATION  OF  DEPENDENT  AND  SYSTEM  PARAMETERS 

The  first  step  to  the  *  term  generation  is  identifying  the  dependent 
parameters  (those  to  be  evaluated  in  the  study).  The  responses  of 
interest  are:  stress  (a),  displacements  (d),  and  acceleration  (a)  due 
to  an  explosive  loading.  The  structure  is  loaded  with  a  point  charge  of 
energy  (En)  with  a  peak  pressure  P0.  The  magnitude  and  characteristics 
of  the  dependent  parameter  on  this  energy  input  and  an  entire  list  of 
system  parameters. 

The  list  of  system  parameters  should  include  all  relevant  material, 
geometric  c.'.d  energy  variables  that  affect  on  the  response  parameters. 
Nielson  (Reference  9)  provided  an  initial  list  of  system  parameters  from 
a  similar  study  involving  a  concrete  burster  slab,  Table  2.  From  the 
standpoint  of  the  GLCM  structure,  additional  parameters  needed  to  be 
included.  Neilson’s  effort  did  not  include  the  material  properties  of 
steel  reinforcement.  Since  the  GLCM  requires  substantial  steel  rein¬ 
forcing,  steel  strength,  cross-sectional  area  and  modulus  were  added  to 
Nielson's  list.  Soil  cohesion,  angle  of  internal  friction,  preconsoli¬ 
dation  pressure,  and  permeability  are  also  added  but  are  easily 
determined  from  Nielson's  work  by  inspection.  Additional  parameters  of 
strain  and  Poisson's  ratio  are  added  but  only  for  clarification  for  they 
are  »  terms  in  themselves.  Table  3  shows  a  complete  list  of  the 
dependent  and  system  parameters  used  in  the  Initial  effort. 
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TABLE  2.  SYSTEM  PARAMETERS  FOR  BLAST  LOADING  (Reference  9) 


Parameter 


Description 


Dimension 


Explosives 


Burster-Si ab 


Energy,  E 
Pressure,  Pp 
Depth  of  Burial,  D 


Thickness,  H 
Mass  Density,  p^ 
Dilatational  Wave  Speed,  Ci 
Poisson's  Ratio, 

Strength, 


ML2T'2 

ML-1T-? 


ML'3 

LT'1 


ML_1T'2 


Soil  (Sand! 


Other 

Parameters 


Mass  Density, 
Dilatational  Wave*  Speed, 
Poisson's  Ratio, 

Strength  Parameter, 


Gravity,  g 

Bucket  Dimensions,  r,  d 
Range 


ML'3 

LT'1 


ML_1T'2 


LT'2 

L 

L 


TABLE  3.  LIST  OF  PROBLEM  PARAMETERS 


DEPENDENT  PARAMETERS 

1.  a  -  STRESS  2.  d  -  DISPLACEMENT  3.  a  -  ACCELERATION 


SYSTEM  PARAMETERS 


I. 

Po 

- 

CHARACTERISTIC  PRESSURE 

14. 

Es 

- 

SOIL  MODULUS 

2. 

I 

- 

IMPULSE 

15. 

c 

- 

SOIL  COHESION 

3. 

En 

- 

ENERGY 

16. 

pc 

- 

PRECONSOLIDATION  PRESSURE 

4. 

R 

- 

RADIUS 

17. 

k 

- 

PERMEABILITY 

5. 

H 

- 

THICKNESS 

18. 

9 

- 

GRAVITY 

A. 

pc 

- 

CONCRETE  MASS  DENSITY 

19. 

T 

- 

TIME 

7. 

Cc 

- 

CONCRETE  P-WAVE  SPEED 

20. 

* 

- 

SOIL  ANGLE  OF  FRICTION 

ft. 

V 

- 

CONCRETE  STRENGTH 

21. 

cst 

- 

STEEL  STRAIN 

9. 

Ec 

- 

CONCRETE  MODULUS 

22. 

Es 

- 

SOIL  STRAIN 

10. 

Fst 

- 

STEEL  STRENGTH 

23. 

t 

c 

- 

CONCRETE  STRAIN 

11. 

Ast 

- 

AREA  OF  STEEL 

24. 

wst 

- 

STEEL  POISSON  RATIO 

12. 

ps 

- 

SOIL  MASS  DENSITY 

25. 

pc 

- 

CONCRETE  POISSON  RATIO 

13. 

cs 

- 

SOIL  P-WAVE  SPEED 

26. 

Est 

- 

STEEL  MODULUS 

8. 

SELECTION  OF  FUNDAMENTAL  UNITS 

The  expression  of  the  dependent  and  system  parameters  in  the  dimen¬ 
sional  matrix  must  be  in  the  form  of  fundamental  units.  These  units  are 
typically  either  mass,  length,  time  (MLT)  or  force,  length,  time 
(FLT).  The  choice  between  e  ch  depends  on  familiarity  of  use  and  ease 
of  »  term  operation.  If  the  w  terms  generated  using  FLT  are  not  easy  to 
satisfy,  changing  to  MLT  may  make  the  scaling  laws  easier  to  satisfy. 
Since  previous  work  used  MLT  (Reference  9)  and  the  fact  that  unit  mass 
does  not  change  in  multiple-gravity  environments,  the  MLT  convention  was 
chosen . 
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C.  SELECTION  OF  REPEATING  PARAMETERS 

To  solve  the  dimensional  matrix,  three  repeating  parameters  must  be 
chosen.  These  parameters  must  be  selected  from  the  list  of  system 
parameters.  The  dependent  parameters  cannot  be  used  without  sacrificing 
the  solution.  The  repeating  parameters  must  be  independent  within 
themselves  to  ensure  independent  solutions.  They  must  contain  all  the 
fundamental  units  at  least  once,  without  being  equal,  in  order  to  ensure 
dimensionless  solutions.  Most  importantly,  they  must  be  easily 
controlled,  experimentally.  This  allows  simple  parameter  modification 
to  meet  scaling  laws  once  derived,  based  on  these  restrictions  impulse 
(I),  height  (H),  and  concrete  modulus  (Ec)  were  chosen  as  the  initial 
repeating  parameters. 
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SECTION  VI 


PI  SET  DETERMINATION 


A.  DIMENSIONAL  MATRIX 

Based  on  the  selection  of  dependent  parameters,  system  parameters 
and  fundamental  units,  a  dimensional  matrix  was  developed.  Table  4. 

Note  that  parameters  23-?6  of  Table  3  are  not  in  the  dimensional 
matrix.  These  are  either  »  terms  in  themselves  or  redundant  due  to 
i nspecti on . 

B.  SOLUTION  MATRIX  AND  it  TERMS 

With  the  dimensional  matrix  and  repeating  parameters  (impulse, 
height,  and  modulus)  identified  in  PISETS,  the  program  generates  the 
solution  matrix.  Table  5.  This  matrix  reduces  to  19  it  terms  as  defined 
by  Equation  (5),  Table  f>.  Note  that  most  of  the  it  terms  are  dependent 
on  the  material  properties  of  the  system  materials.  Contingent  on  the 
use  of  the  same  material  in  the  model  ana  prototype, the  resulting 
scaling  laws  fall  out  quite  easily. 

C.  EVALUATION  OF  it  TERMS  AND  SCALIN'!  LAWS 

The  simplest  way  to  organize  an  evaluation  process  Is  to  group  the 
generated  »  terms  into  functional  groups  relating  their  driving 
influence  on  the  system.  The  »  terms  are  thus  separated  into  geometric, 
strength,  explosive  and  impulse  groups.  In  this  way  the  impact  of  a 
violated  »  term  on  the  problem  can  be  readily  identified.  For  Instance, 
if  a  *  term  defining  a  strength  relationship  requires  the  modulus  of 
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Repeating  Variables 


TABLE  6.  INITIAL  GENERATED  n  TERMS 
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steel  to  be  twice  that  of  the  prototype,  its  effect  cn  other  strength 
relationship  is  relatively  identifiable.  Table  7  indicates  the 
appropriate  grouping. 


TABLE  7:  it  TERM  GROUPS 


GROUP 

w  TERMS 

Geometric 

*2’  "6*  *11 

Strength 

*9’  *10’  *14 *  *15’  *16 

Expl osi ve 

*!»  n4>  ^5 

Impulse 

*3’  V  V  *12 

13’  17’  18’  *19 

Once  grouped,  the  actual  scaling  laws  are  formulated  by  equating 
model  it  tern  to  prototype  »  tern  (Equation(4)),  using  the  length  scale 
factor  n  (Equations  (27)  and  (20)), and  cancelling  out  like  material 
properties.  The  geometric  group  results  in  the  following,  where 
subscripts  m  and  p  are  model  and  prototype,  respectively: 
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This  n  set  developed  with  impulse,  height,  and  modulus  of  elasticity 
satisfies  all  three  requirements.  Equations  (46),  (47),  (48),  (53), 

(55),  (56),  (57),  (63),  and  (64)  are  all  consistent  with  Table  1. 

Equation  (49),  (50),  (51),  (52),  (54),  (58),  (59),  (60),  (61),  and  (62) 
are  consistent  with  the  assumption  that  material  properties  between  model  anc 
prototype  are  the  same.  It  appears  that  all  w  terms  are  physically 
consistent  as  well;  however,  a  re-evaluation  of  the  system  parameters 
indicates  a  more  applicable  w  set  may  be  available. 

The  present  list  of  system  parameters  in  the  dimensional  matrix 
consists  of  19  quantities.  Based  on  GLCM  site  definition  and  redundancy 
within  the  list,  it  is  author's  contention  that  only  17  are  needed 
without  losing  clarity  or  applicability.  The  two  parameters  dropped  are 
impulse  and  permeability.  Impulse  is  a  function  of  pressure,  area,  and 
time.  Since  all  three  are  defined  in  the  system  parameter  list,  impulse 
is  redundant.  With  regards  to  permeability,  it  is  only  reievanc  when 
water  or  some  pore  fluid  is  present.  Since  GLCM  sites  do  not  have 
swallow  water  tabl es,  permeahi 1 i ty  is  not  a  concern.  Based  on  these 
considerations  a  new  *  set  should  be  determined  from  the  reduced  system 
parameter  list. 

D.  ALTERNATIVE  it  SET 

By  deleting  impulse  the  previous  choice  of  repeating  parameters 
cannot  be  maintained.  A  new  selection  is  required.  This  leads  to  a 
trial  and  error  process, whereby  the  following  methods  can  be  used  to 
find  a  reasonable  it  set. 

1.  Combine  the  existing  tt  termu  by  multiplication  or  division. 

2.  Select  new  repeating  variables 
or 

3.  Change  the  system  of  fundamental  unit  (MLT  or  FLT) 
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With  the  existence  of  PISETS  the  simplest  method  is  to  select 
another  sec  of  repeating  variables  and  simply  re-evaluate  the 
new  it  set.  With  each  new  selection,  PISETS  will  check  for  independence 
and  calculate  a  solution  matrix.  The  operator  must  ensure  that  the 
otner  repeating  variable  requirements  listed  in  Section  II  are  met. 
Several  new  it  sets  were  considered.  The  preferred  it  set  is  discussed  in 
Section  VII. 


SECTION  VII 


FINAL  w  SET  DETERMINATION 


A.  DIMENSIONAL  MATRIX 

The  preferred  u  set  for  the  GLCM  shelter  resulted  from  the  selection 
of  concrete  modulus  of  elasticity  (Ec),  concrete  mass  density  (oc)  and 
fundamental  height  (H)  as  the  repeating  variables.  Table  8  indicates 
the  dimensional  matrix  used  in  the  final  solution. 

Note  that  impulse  and  permeability  are  no  longer  in  the  list  of 
system  parameters.  Impulse  was  redundant,  as  discussed  in  Chapter  VI, 
and  was  therefore  deleted.  Permeability  was  deleted,  due  to  the 
unlikelihood  of  a  high  water  table  at  potential  GCLM  shelter  sites. 

This  is  often  a  favorable  deletion  because  of  the  difficulty  of 
realistically  evaluating  permeability  in  a  centrifuge.  The  use  of  the 
same  soil  in  model  and  prototype  does  not  mean  perm  abilities  will  be 
the  same.  Although  the  pore  passages  will  be  the  s..me,  the  driving 
forces  are  not;  therefore  neither  is  permeability.  This  will  be  further 
discussed  later  in  this  section. 

8.  SOLUTION  MATRIX  AND  »  TERMS 

Inputting  the  dimensional  matrix  and  identifying  Ec,  pc,  and  H  as 
the  repeating  variables  in  PISETS,  the  solution  matrix  Is  generated, 
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Repeating  Variables 


Note  that  in  comparison  to  Table  5, the  diagonal  number  seauence 
across  the  matrix  is  2  and  not  1.  This  indicates  that  in  generation  of 
the  homogeneous  algebraic  equation  for  the  repeating  variables  similar 
to  Equations  (15),  (16),  and  (17),  fractions  resulted  in  the 
exponents.  P I  SETS  does  not  operate  with  fractional  exponents  and 
therefore,  multiplies  the  exponent  appropriately  in  order  to  get  an 
integer.  In  this  case  each  term  was  squared.  When  the  it  terms  are 
actually  written, the  square  root  can  be  taken  to  bring  them  back  to  the 
lowest  set  of  integer  exponents.  For  instance 


IT 


1 


which  reduces  to 


(65) 


=  jr-  (66) 

while  remaining  dimensionless. 

The  matrix  results  in  17  w  terms,  two  less  than  the  initial  trial 
matrix  in  Table  10.  This  reduction  is  consistent  with  Equation  (5)  and 
the  reduction  in  system  parameters. 

C.  EVALUATION  OF  *  tpRHS 

As  with  trial  it  sets,  this  *  set  is  grouped  for  evaluation.  It  is 
divided  into  geometric,  strength,  explosive  and  mass  density  groups, 
shown  in  Table  11. 
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TABLE  10.  SOLUTION  »  TERMS 
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TABLE  11:  SOLUTION  ir  TERM  GROUPS 


GROUP 

it  TERMS 

Geometri c 

V  V  *10 

Strength 

*8’  V  *13*  *14 ’  *15 

Explosive 

V  *4’  *5 

Mass  Oensity 

*3 ’  *7’  *11’  *12’  *16’  *17 

The  groups  are  reduced  to  the  scaling  laws  in  the  same  way  as  the 
initial  set  (Section  VI).  Tne  geometric  group  resulted  in  the  following 
laws: 
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Since  the  strength  group  contains  only  material  property  ratios,  using 
the  same  materials  in  model  as  prototype  results  in  a  one  to  one 


rel aticnship.  Therefore,  the  strength  group  dictates  the  following 
scaling  laws: 
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The  explosive  group  reduces  to  three  scaling  laws. 


(70) 

(71) 

(72) 

(73) 

(74) 


o 

*1  =  F" 
c 


c 

_  r  m\ 

'tjT 

P 


a 

P 


.* .  a  «  a 
m  p 


(75) 


46 


4  '  r~ 

c 


p„  p 

or  0 

P  _  n 


Po„  ■  <r=>  \ 

m  c  r 


•••  p0  -Po  (76) 

m  p 


En 

p  _  m 


K  p  £c  H'n  Ec 
p  m 


H3  E 
me  r 
c  ,  f  mi  E 

nm  H3”  E  np 

P  CP 


*'En  "-4-  (77) 


m  n 


The  mass  density  group  reduces  as  folic 
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Cnee  determined,  the  scaling  laws  are  checked  for  consistency  with  the 
three  scaling  requirements  listed  in  section  6.3 

Equation'  (67),  (68),  (69),  (75),  (76),  (77),  (78),  (79),  (82)  and 
(83)  are  all  consistent  with  Table  1.  They  indicate  that  the  standard 
scaling  laws  are  not  violated  and,  therefore,  satisfy  the  first  condition 


/or  a  workable  *  set.  Equations  (70),  (71),  (72),  (73),  (74),  (79),  and 
(81)  are  consistent  with  the  use  of  the  sane  material  properties  between 
model  and  prototype  and  satisfy  the  second  condition  for  a  viable  *  set. 
The  last  condition  is  more  difficult  to  determine.  Are  the  scaling  laws 
physically  possible  to  satisfy  and  maintain  within  themselves? 

The  answer  to  this  question  leads  to  several  inconsistencies  and 
conflicts  between  assumed  properties  and  derived  scaling  laws.  The 
fo1  lowing  discussion  does  not  invalidate  this  ir  set  but  qualifies  its 
use.  Note  at  the  outset  other  researchers  (References  8,14)  have 
identified  these  problems  and  have  either  suggested  solutions  or 
dismissed  the  problems  away. 

The  major  problem  stems  from  the  internal  conflict  between 
satisfying  the  geometric  scaling  laws  while  maintaining  the  material 
properties.  Just  as  the  external  length  of  the  GLCM  structure  must  be 
reduced  by  n,  according  to  Equation  68,  the  internal  dimensions  of  the 
.concrete  and  soil  must  be  as  well.  For  soil,  a  reduced  grain  size 
drastically  changes  its  material  properties.  For  concrete,  similar 
material  properties  pertubations  occur  (Reference  15).  For  concrete,  an 
answer  is  provided  by  Sabnis  and  White  (Reference  14),  who  show  that 
gypsum  mortar  can  be  used  instead  of  concrete  on  small-scale  models. 
Their  work  indicates  that  prototype  material  properties  are  duplicated 
while  satisfying  geometric  scale.  Note  that  the  gypsum  particle  may  not 
exactly  meet  the  scaled  particle  requirements  but  icr7  meter  size 
particles  do  not  cause  perturbations  and  therefore  can  be  used. 

For  soil,  the  solution  is  not  as  simple.  Soil  is  a  two-phase  material 
consisting  of  solids  and  fluids.  The  various  interactions  of  these 
phases  make  it  impossible  to  change  one  property  without  affecting 


50 


others.  If  the  particle  size  is  reduced, all  of  the  m  terms  relating 
soil  modulus,  mass  density,  and  dilational  wave  speed  change.  It  is 
appropriate  to  avoid  these  alterations  by  ignoring  the  geometric 
reduction  for  soil  particles. 

By  using  the  same  soil  in  model  as  in  prototype  consistent  material 
properties  can  be  assumed.  Any  question  as  to  adverse  effects  on  soil 
structure  interaction  due  to  this  procedure  was  addressed  by  Morris 
(Reference  8).  He  indicates  that  if  the  model  structure  is  sufficiently 
large  with  respect  to  soil  particle  size,  the  material  properties  can  be 
maintained  without  significant  effect  on  structure  soil  interaction. 

This  sufficient  size  differential  will  change  as  the  model  changes.  To 
determine  the  limiting  differential,  a  modeling  of  models  procedure 
should  be  conducted.  This  can  be  done  by  testing  a  generic  model  at 
different  gravities  and  checking  for  inconsistencies.  If  this  is  done 
and  the  scale  chosen  accordingly,  the  same  soil  can  be  used  without 
concern. 

The  secondary  problem  of  the  *  set  was  alluded  to  in  Section  VII  A. 

It  is  common  in  all  centrifuge  testing,  using  the  same  soil  as 
prototype.  Ic  is  often  ignored  as  in  the  generation  of  this  it  set,  but 
it  needs  to  be  clarified.  The  problem  is  permeability.  On  exposure  to 
multiple-gravity  forces,  water  will  pass  through  soil  pore  passages 
faster  than  at  one  gravity.  The  hydrodynamic  scaling  law  (Equation  (44)) 
indicates  this  fact.  When  stresses  build  up  in  a  structure/soil  system 
in  which  water  is  present,  the  stresses  in  the  structure  will  dissipate 
at  n  times  faster  than  the  prototype  but  the  pore  water  pressure  will 
dissipate  at  n^  time  faster.  This  negates  similitude  between  model  and 
prototype.  It  may  be  possible  to  introduce  a  more  viscous  fluid  into 
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the  problem  to  slow  the  dissipation,  but  this  would  require  an 
evaluation  of  the  corresponding  changes  in  soil  shear  strength.  It  is 
sufficient  to  be  aware  of  this  problem  and  avoid  it  by  keepiny  water 
tables  out  of  the  system. 

Based  on  these  assumptions  and  clari fications ,  the  third  requirement 
is  satisfied.  The  scaling  laws  are  physically  consistent  and  possible 
within  themselves.  From  the  standpoint  of  internal  consistency  of  the 
scaling  laws  and  similitude,  the  stage  is  set  for  the  model  testing  of 
the  (1LCM  shelter. 


SECTION  VIII 


PRACTICAL  LIMITATIONS 


A.  LIMITATIONS 

With  the  basic  scaling  laws  established,  the  scale  to  which  the  GLCM 
model  can  actually  be  tested  must  be  ascertained.  The  applicable  scale 
to  be  used  must  be  evaluated  outside  of  the  scaling  laws  and  theory  of 
multiple-gravity  testing.  The  physical  limitations  of  equipment, 
structure  constructability,  soil  sample  preparation,  explosive,  and 
instrumentation  play  a  vital  role  in  the  determination  of  scale.  Each 
of  these  considerations  may  or  may  not  induce  uncertainties  into 
testing,  depending  on  the  scale  and  complexity  of  the  problem.  The  GLCM 
shelter  testing  program  is  susceptible  to  limitations  by  all  five 
considerations. 

1.  Equipment 

The  primary  factor  dictating  the  selected  scale  or  size  of  the 
model  is  the  availability  of  centrifuges  built  to  test  dynamic  behavior. 
Presently,  within  the  United  States,  eight  existing  centrifuges  and  two 
potential  centrifuges  are  suitable  for  a  GLCM  testing  program.  Table  12 
indicates  their  location  and  respective  capabilities.  The  recommended 
model  tieights  and  lengths  are  based  on  Avgherinos'  and  Schofield's 
(Reference  16)  work  and  Bassett's  (Reference  17)  work.  Because  of  the 
radical  nature  of  centrifugal  accelerations  applied  to  a  model  within  a 
centrifuge,  the  gravitational  forces  on  the  model  will  vary  across  its 
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base  and  height,  unlike  that  of  the  prototype.  To  limit  this  effect, 
Avgherinos  and  Scho'ieM  (Reference  16)  suggest  that  the  height  of  the 
model  be  limited  to  1/10  of  the  radial  arm  of  the  centifuge.  Bassett 
(Reference  17)  suggests  the  length  be  no  longer  than  the  cord  of  a  15° 
section  with  a  radius  equal  to  the  centrifuge's  arm.  The  values 
provided  are  theoretical  and  do  not  reflect  the  actual  physical 
requirements  of  the  machine  itself. 

Aside  from  the  acceleration  variation,  two  additional  equipment 
concerns  affect  the  choice  of  scale.  The  scales  suggested  in  Table  12 
do  not  consider  the  weight  limit  associated  with  the  maximum  centrifugal 
accel erations .  As  the  payload  weight  increases,  the  maximum  available 
gravities  decrease.  If  the  eventual  GLCM  model  weighs  more  than  the 
specific  weight  limit,  the  corresponding  decrease  in  acceleration  must 
be  calculated  into  the  required  scale.  Secondly,  there  is  no  allowance 
for  edge  effects.  As  the  explosive  energy  strikes  the  base  or  sides  of 
the  spinning  platform,  it  reflects  back  into  the  model,  disrupting  the 
model  response.  Space  should  be  allowed  between  the  model  and  these 
surfaces  to  minimize  these  effects.  Consequently,  the  actual  scale  may 
have  to  be  smaller  than  indicated  in  Table  12.  The  actual  space  needed 
will  depend  on  the  energy  absorbability  of  the  spinning  platform 
surfaces,  tne  energy-transmi ssion  ability  of  the  soil  and  the  length  of 
response-time  histories  required.  A  recommended  spacing  will  have  to  be 
evaluated  experimental ly  for  each  centrifuge  used.  An  experiment  using 
just  a  soil  mass,  loaded  by  a  small  explosive  surface  charge,  should 
suffice  in  providing  the  information. 
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2.  Constructability  of  Model  Structures 

To  satisfy  the  scale  needed  to  actually  place  a  model  in  a 
centrifuge,  the  structure  must  be  constructabl e.  As  previously 
discussed,  gypsum  mortar  has  been  suggested  oy  Sabnis  and  White 
("Reference  14)  for  the  structural  material.  The  extent  to  which  it  can 
be  poureo  to  a  1/100  and  smaller  scale  is  not  known.  For  example,  if 
the  New  Mexico  Engineering  Research  Institute's  (NMERI)  centrifuge  was 
used,  the  interior  wall  thickness  of  the  GLCM  shelter  model  would  have 
to  be  at  least  0.2  inches  thick.  When  considering  that  steel 
reinforcement  has  to  be  placed  within  that  wall,  the  difficulty  of  the 
problem  becomes  apparent.  Although  there  is  some  research  in  scale 
model  reinforcement  (Reference  18)  this  effort  will  require  an 
independent  study.  The  walls  are  too  small  with  respect  to  the  amount 
of  reinforcement  required. 

3.  Model  Soil  Preparation 

The  prototype  soil  material  is  primarily  granular.  Tivs 
significantly  reduces  the  complexity  of  the  model  soil  preparation. 

Size  effects  present  no  significant  problem  if  Morris'  (Reference  8) 
guidance  (Section  VII  C)  is  adhered  to.  Consolidation  effects, 
primarily  associated  with  clay  are  essentially  nonexistent.  The 
previously  discussed  inconsistency  with  permeability  is  avoided, 
assuming  water  is  excluded.  Reconstruction  of  stability  within  the  soil 
from  chrmical  processes  is  of  no  concern  since  the  prototype  requires 
disturbed  soil.  Therefore,  it  is  sufficient  for  testing  to  place  the 
soil  material  at  the  requireo  density  and  allow  the  appropriate  stresses 
to  develop  during  the  actual  testing  of  the  model  (Reference  19). 
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4.  Explosives 

8y  virtue  of  the  scaling  lav/s  established,  energy  scales  at  l/nJ 
that  of  the  prototype  energy.  This  requires  very  small  energy  sources 
for  excitation  of  the  model  structure.  For  instance,  the  NMERI's 
centrifuge  would  need,  at  a  minimum,  a  charge  of  .0007  pounds  of  TNT  to 
match  a  2,000-pound  prototype  threat.  This  is  1/100  of  an  ounce  of  TNT 
or  equivalent  explosive.  For  conventional  threats  of  500-pound  to 
2,000-pound  bombs  this  could  be  very  difficult.  For  nuclear  threats  the 
downsizing  of  the  energy  source  would  be  less  difficult  due  to  several 
orders  of  magnitude  increase  in  prototype  source.  Regardless,  the 
equivalent  explosive  is  quite  small  and  presses  the  techrical  limits  of 
commercial  explosives. 

Nielson  (Reference  9)  has  proposed  that  cyclothimethlenetrini - 
tromine  (RDX)  be  used  for  the  above  expressed  purpose.  Unfortunately, 
his  efforts  do  not  go  down  to  the  scales  needed  for  the  GLCM,  Table  13. 
However,  with  a  modification  to  a  Reynold  RP-83  detonator,  he  has 
recorded  some  reasonable  data  about  the  energy  source  itself.  It  may  be 
possible  to  reduce  this  detonator  even  more  for  GLCM  application. 

It  is  generally  appropriate  to  apply  scaling  laws  to  the 
explosive  shape,  as  well  to  its  energy  output.  Nielson  suggests  that 
the  majority  of  the  randomness  in  his  data  is  attributable  to  explosive 
geometric  variation  and  violation  of  the  geometric  scaling  law  of  1/n. 
But  since  the  GLCM's  threat  is  projectile  in  nature  and  it  is  surrounded 
by  a  soil  overburden,  the  geometric  violations  may  not  be  a  large 
consequence.  If  the  peak  pressure  and  duration  of  the  explosive  can  be 
appropriately  duplicated  by  the  time  it  strikes  the  structure,  the  shape 
requirements  may  possibly  be  ignored.  At  present,  it  Is  unclear  whether 
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this  simplification  or  explosive  can  he  used  on  a  model  GLCM  shelter. 
If  the  proposed  model  testing  is  to  be  limited  to  the  available 
centrifuges  an  explosive  evaluation  and  development  program  will  be 
required.  Considering  several  potential  threats  to  the  GLCM  shelter, 
this  could  be  a  costly  development. 


TABLE  13.  THEORETICAL  MODEL  EXPLOSIVE  SIMULATION  WEIGHTS 

(Reference  9) 


Threat 

Designation 

(lbs) 

Centri fuge 

Envi ronment 

(Gravities ) 

20  g 

40  g 

60  g 

80  g 

100  g 

Weight 

of  RDX  in 

grams 

250 

12.35 

1.54 

0.46 

0.19 

0.10 

500 

25.57 

3.20 

0.95 

0.40 

0.20 

1000 

53.43 

6.68 

1.98 

0.83 

0.43 

2000 

106.95 

13.37 

3.98 

1.67 

0.86 

5.  Instrumentation 

Like  other  limitations  on  physical  equipment  of  a  potential 
testing  system,  instrumentation  is  a  consideration.  With  small -scaled 
structures,  the  placement  of  instrumentation  cannot  inhibit  the 
structure's  response.  They  must  be  small  enough  and  be  placed  in  such 
manner  and  number  as  to  prevent  perturbations  in  response.  Fortunately 
the  technology  exists  to  meet  the  GLCM  testing  needs.  Sabnls  et  al. 
(Reference  13)  indicate  that  electrical  resistance  strain  gages  exist  as 
small  as  1/64  of  an  Inch  thick.  Ryan’s  (Reference  20)  work  on  piles 
used  such  gages  with  reasonable  success. 
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For  stresses,  two  possibilities  exist,  direct  or  indirect 
measurement.  Nielson  (Reference  9)  suggests  Dynasen  carbon  foil 
pressure  sensors  as  a  direct  measure.  Others  suggest  indirect  measures 
using  Quartz  transducers  to  measure  velocity  and  then  calculating  stress 
by  Equation  (83). 

a  =  pcv  (84) 

where  p  =  density,  c  =  speed  of  light  and  v  =  measured  velocity 
(Reference  21). 

Other  dynamic  quantities  such  as  acceleration  and  displacement 
can  be  measured  using  variations  on  the  velocity  quartz  transducer. 
Complete  specifications  can  be  obtained  through  manufacturers  listed  in 
Table  14.  Note:  Although  the  gages  exist  for  the  quantities  to  be 
measured,  caution  should  be  used  in  the  development  of  the 
instrumentation  plan  concerning  space  limitations  for  gage  connections 
to  external  recording  equipment. 
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TABLE  14.  PARTIAL  LIST  OF  INSTRUMENTATION  MANUFACTURERS 
(References  9,  13) 


Manufacturer 


Acurex  Corporation 


BLH  Electronics 


Dynasen,  Inc. 


Mi cromeasurements 
Div.  of  Vishay 
Inter-Technology, 
Inc. 

PCB  Peizotrorics, 


Strainsert  Company 


Address  Type  of  Accessories 


485  Clyde  Avenue 
Mountain  View,  C A 
94042 

Tel:  (415)  964-3200 

42  Fourth  Avenue 
Waltham,  Mass.  02154 
Tel:  (617)  890-6700 

20  Dean  Arnold  Place 
Goleta,  CA  93017 

P.0.  Box  27777 
Raleigh  NC  27611 
Tel:  (919)  365-3800 


P.C.  Box  33 
Buffalo  NY  14225 
Tel:  (716)  684-0001 

Union  Hill  Industrial 
Park 

West  Conshocken, 

PA  19428 


Data  acquisition  system 
for  strains 


All  types  of  strain  gages 
and  related  accessories 
for  strain  reading 

Pressure  Sonsors 


All  types  of  strain  gages 
and  related  accessories 
for  strain  reading 


Quartz  Transducers  for 
dynamic  measurements 


Flat  load  cells  of 
various  capacities 
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SECTION  IX 


RECOMMENDATIONS 


A.  SCALING  LAWS 

The  technical  and  financial  benefits  of  the  centrifuge  are  set.  The 
extent  to  which  these  benefits  can  be  realizes  're  not  and  are  dependent 
on  achieving  similitude  between  the  model  and  prototype.  Using  the 
scaling  laws  in  Table  15,  this  dependency  can  be  satisfied  and  a  viable 
testing  program  can  result.  If  adhered  to,  these  scaling  laws  will 
creditably  relate  th_-  GLCM  model  shelter  response,  point  for  point,  to 
the  prototype.  These  laws  have  been  chosen  to  meet  the  standard  scaling 
laws  associated  with  common  material  properties,  internal  consistency 
and  centrifuge  capabilies;  therefore,  the  extent  to  wKich  the  laws  are 
satisfied  will  depend  on  the  skill  of  researcher  and  the  centrifuge 
chosen. 

B.  MATERIALS 

Based  on  the  available  data  and  the  considerable  simplification  of 
the  modelling  process,  the  materials  used  in  model  end  prototype  should 
have  the  same  material  properties. 

Gypsum  mortar  is  suggested  to  be  used  as  a  concrete  substitute 
because  it  exhibits  the  appropriate  material  properties  while  avoiding 
aggregate  scaling.  The  reinforcement  should  be  of  the  same  quality  and 
strength  as  the  prototype.  Whereas  soil  qiality  and  strength  should  be 
exactly  the  same  as  prototype  to  avoid  mater. ’1  property 
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TABLE  15:  RECOMMENDED  SCALING  LAWS 


PARAMETER 


SYMBOL 


SCALING  LAW 


1.  Stress 

a 

0,  = 

m  p 

d 

2.  Displacement 

d 

d 

m  n 

3.  Acceleration 

a 

am  =  n  ap 

4.  Velocity 

V 

vm  =  vp 

5.  Explosive  Pressure 

po 

P  =  P 
°m  °p 

En 

6.  Explosive  Energy 

En 

R 

En  “ 
m  nJ 

7.  Radius 

R 

R_  =  -E. 

m  n 

H 

8.  Thickness 

H 

Hm  =  -E 

m  n 

9.  Material  Density 

p 

p  = 

m  p 

10.  Material  Modulus 

E 

Em  ~  Ep 

11.  Material  Strength 

F 

Fm  =  Fp 

12.  Material  Dilitional 

C 

wave  speed 

A 

13.  Area 

A 

n 

V 

14.  Volume 

V 

M 

v-'i 

15.  Mass 

M 

n 

16.  Material  Strain 

e 

« * 

c  »  e 

t 

m  p 

17.  Dynamic  Time 

t 

t  - -E. 
n 

18.  Material  Poisson's  Ratio 

u 

u  »  u 

m  p 

19.  Soil  Cohesion 

c 

cm  *  cp 

20.  Soil  Preconsolidation 

Pc 

pc  “  pc 

Pressure 

ffd 

cm  cp 

21.  Force 

ff 

ff  »-/ 

m  rr 

22.  Soil  Angle  of  Internal 

♦ 

♦_  “  ♦_ 

Friction 

m  p 
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inconsistencies.  Any  deviation  from  this  recommendation  will  cause 
numerous  inconsistencies  in  the  scaling  laws  and  several  physical 
impossibilities  in  material  behavior  and  experimental  capabilities.  The 
sources  of  uncertainty  due  to  equipment  and  model  construction  are 
numerous  enough  without  introducing  material  property  variations.  To 
avoid  the  inconsistency  of  the  rate  of  structure  and  pore  pressure 
stress  dissipation  due  to  changes  in  permeability,  water  should  be  left 
out  of  the  problem.  Assuming  a  low  water  table  at  potential  GLCM  sites, 
this  restriction  is  consistent  with  the  prototype. 

C.  CENTRIFUGE 

The  centrifuge  chosen  for  the  GLCM  shelter  model  testing  controls 
the  scale  to  be  chosen  for  testing.  It  is  beneficial  to  the  success  of 
the  program  to  choose  the  largest  centrifuge  available.  As  the  scale 
gets  larger  the  complexity  and  likelihood  of  experimental  error  reduces 
due  to  material  and  equipment  limitations.  The  author  believes  that  the 
complexity  of  the  GLCM  requires  the  use  of  the  national  centrifuge 
research  facility  at  Ames  Laboratory.  This  could  possibly  limit  the 
scale  reduction  to  1/38.  The  use  of  any  smaller  facility  will  require 
extensive  simpl icatior.  of  the  model  design  due  to  model  construction 
problems. 

D.  E/PLOSIVES 

The  availability  of  controllable  explosives  of  the  sizes  needed  for 
this  model  study  is  limited  at  best.  Nielson's  (Reference  9)  work  using 
RDA  (cyclothimethlenetrinitromine)  is  encouraging  but  not  conclusive. 

If  large  centrifuges  similar  to  Ames  facility  are  not  available,  a 
research  program  to  develop  an  energy  source  small  enough  to  be 
consistent  with  GLCM  threats  should  be  conducted. 


E.  INSTRUMENTATION 

Instrumentation  to  measure  stress,  strain,  and  dynamic  quantities  is 
available  commercially.  The  exact  gages  will  have  to  be  determined  once 
the  centrifuge  and  scale  are  determined.  The  number  of  gages  per  model 
will  depend  on  the  recording  capability  of  the  facility  and  wiring 
capabilities  of  the  testing  technician.  The  selection  of  gages  shoulo 
come  from  manufacturers  listed  in  Table  14. 

F.  ASPECTS  TO  BE  EVALUATED 

Before  actually  testing  the  potential  GLCM  shelter  model,  several 
aspects  of  the  testing  setup  should  be  evaluated.  To  adequately 
establish  the  testing  program,  a  series  of  initial  tests  must  be  done  to 
answer  the  following  questions: 

1.  To  what  extent  must  measures  be  taken  to  limit  edge  effects? 

2.  To  what  scale  must  we  limit  the  testing  to  avoid  excessive 
perturbations  on  the  test  data  from  inappropriate  soil  structure 
interaction  caused  by  unsealed  soil  grain  size? 

3.  What  is  the  optimum  explosive  charge  shape  to  use? 

4.  To  what  extent  do  strain  rate  effects  disturb  this  particular 
structure? 

5.  To  what  detail  can  we  build  the  GLCM  shelter  model,  based  on 
physical  construction  limitations? 

G.  MODELLING  OF  MODELS 

After  the  questions  have  been  answered  and  scale  set  and  testing 
procedure  have  been  determined,  the  setup  should  be  exposed  to  a 
modelling  of  models.  This  will  not  require  a  model  for  each  gravity 
environment  selected.  Generic  setups  can  be  used  to  check  consistency 
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of  developed  procedures.  A  series  of  environments  such  as  20,  40,  60, 
80,  and  100  gravities  should  indicate  any  inconsistency  of  the  testing 
procedure  due  to  the  multiple-gravity  loadings. 

H.  CONSTITUTIVE  MODELLING 

Because  of  scaling  limitations  (GLCM  shelter  size  and  capabilities 
of  existing  U.S.  centrifuges,  and  uncertainties  concerning  edge  effects, 
strain  rate,  instrumentation,  and  explosive  shape  and  composition) 
multiple-gravity  loadings  of  GLCM  shelter  components  are  recommended. 
These  test  results  should  be  analyzed.using  available  numeric  computer 
codes  to  ascertain  compatibility  between  computer  constitutive  modelling 
and  centrifugal  modelling.  If  acceptable  computer  predictions  of 
observed  centrifugal  model  events  can  be  made,  then  centrifugal 
modelling  of  the  entire  GLCM  shelter  may  not  be  required. 
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A.  INTRODUCTION 

The  PISETS  workspace  contains  functions  which  are  used  to  generate 
all  possible  sets  of  nondimensional  parameters  for  a  given  problem  in 
dimensional  analysis.  If  the  generation  of  all  sets  is  not  desirable,  a 
single  set  or  group  of  sets  may  be  specified. 

The  user  must  supply  as  input  the  total  number  of  basic  dimensions 
used  and  the  total  number  of  dimensioned  parameters  used.  The  user  must 
also  develop  and  input  the  dimensional  matrix  that  specifies  which  basic 
dimensions  appear  in  each  dimensioned  parameter. 

Output  will  be  in  the  form  of  a  matrix  of  exponents  which  will 
indicate  to  what  power  each  dimensioned  parameter  must  be  raised,  prior 
to  multiplication,  to  form  each  nondimensional  parameter.  An  exponent 
matrix  will  be  generated  for  each  set  of  nondimensional  parameters. 

PISETi  was  implemented  on  the  Commodore  SuperPET  microcomputer  using 
Waterloo  Computing  System  microAPL  version  1.0. 

PISETS  was  Implemented  on  the  Northeast  Regional  Data  Center  (NERDC) 
system  through  the  Virtual  Machine/System  Product  Conversational  Monitor 
System  (VM/SP  CMS)  using  VS  APL  (APL*PLUS/1140) . 

Commands  to  be  entered  by  the  user  will  be  indicated  1  t  boldface 
throughout  this  manual. 


71 


B.  THE  DIMENSIONAL  MATRIX 


A  dimensional  matrix  for  tiie  problem  should  be  developed  prior  to 
the  use  of  P I  SETS . 

Each  column  of  the  dimensional  matrix  corresponds  to  one  of  the 
dimensioned  parameters;  each  row  corresponds  to  one  of  the  basic 
dimensions.  The  values  in  the  matrix  represent  the  power  to  which  each 
basic  dimension  appears  in  each  dimensioned  parameter. 

As  an  example,  consider  a  problem  for  which  three  basic  dimensions 
are  chosen;  Mass  (M),  Length  (L)  and  Time  IT).  Assume  that  the  problem 
(a  spri ng-mass -dashpot  system)  is  fully  described  by  the  seven 
dimensioned  parameters,  with  their  component  basic  dimensions,  shown 


below : 

y,  displacement,  L 

y0,  initial  displacement,  L 

n ,  mass,  M 

c,  damping  coefficient,  MT~* 

k,  spring  constant,  MT~^ 

vQ,  initial  velocity,  LT~* 

t,  time,  T 


The  dimensioned  parameters  are  formed  from  powers  of  the  bas’c 
dimensions.  These  exponents  are  entered  into  the  dimensional  matrix,  as 
shown  below: 


1 

2 

3 

4 

5 

6 

7 

y 

m 

c 

k 

vo 

*0 

t 

L 

l 

0 

0 

0 

1 

0 

T 

0 

0 

-I 

-2 

-1 

0 

1 

M 

0 

1 

1 

1 

0 

0 

0 
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It  is  important  to  record  the  column  number  that  corresponds  to  each 
dimensioned  parameter.  When  the  exponent  matrix  is  output  each  columr 
of  that  matrix  will  correspond  to  the  same  dimensioned  parameter  that 
the  same  column  of  the  dimensional  matrix  corresponds  to. 

C.  HOW  TO  START  ANO  STOP 

1.  The  Commodore  Super PET  Microcomputer 

Turn  on  both  the  SuperPET  and  the  disk  drives  and  set  the 
SuperPET's  PROG: 6502: 6809  switch  to  6809  and  PROG:R/W:READ  switch  to 
R/W.  The  screen  will  then  display  the  Waterloo  microsystem  menu. 

Put  the  Waterloo  System  diskette  in  disk  drive  1  and  the  PISETS 
diskette  in  disk  drive  0. 

From  the  Waterloo  microSystem  menu,  select  API  by  entering  the 

command: 

apl 

After  a  minute  or  so,  the  screen  will  display  a  message  similar 
to  the  following: 

WATERLOO  MICRO  APL  VERSION  1.0  81/09/01 

COPYRIGHT  1981  BY  WATERLOO  COMPUTING  SYSTEMS  LIMITED 

CLEAR  WS 

Load  the  PISET  workspace  by  entering  the  command: 

>  LOAD  PISETS 

The  screen  will  respond  with  the  date  and  time  that  the 
workspace  was  last  saved.  The  PISET  workspace  Is  now  ready  for  use. 

To  stop,  enter  the  command: 

>  OFF 

and  remove  the  disks  from  the  disk  drives. 


2)  VS  APL  under  NERDC 


An  example  of  the  start  of  a  VS  APL  terminal  session  is  shown  on 

page  82. 

A  *-erminal  with  the  APL  character  set  is  required  when  using  VS 
APL.  If  a  deckwriter  is  used  you  will  have  hardcopy  of  your  output.  To 
get  the  terminal's  attention  you  must  push  the  BREAK  key  and  then 
carriage  return.  Often,  this  must  be  done  repeatedly. 

The  terminal  will  eventually  respond  with  the  prompt: 

enter  class 
Respond  by  entering: 
nerdc 

If  NERDC  is  available  the  terminal  will  print  a  system  menu  (do 
not  select  1  for  APL,  it  is  a  different  system)  from  which  vm/cms  may  be 
selected  by  entering: 

v 

The  terminal  will  print: 

vm/370  c  nl i ne 

Push  carriage  return  and  enter  the  logon  command  as  follows: 

logon  xx/xxxxx,yy 

where  xxxxxxxx  represents  your  access  number  and  yy  is  your  sequence 
number.  The  terminal  will  respond  with: 

ENTER  PASSWORD: 

At  this  point  you  must  type  in  your  password. 

Next  a  block  of  times  and  messages  will  be  printed  ending  with: 
NERDC  VM/SP3  CMS 

Indicating  that  you  are  now  in  the  CMS  system  but  not  yet  in  VS  APL. 
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To  access  APL  enter  the  command: 


apl 

A  block  of  statistics  will  be  printed  followed  by: 

APL*PLUS/li40 

If  a  C0N1INUE  workspace  was  saved  previously,  it  will  now  be 
loaded  by  the  system  and  the  date  and  time  that  this  workspace  was  saved 
will  be  printed.  The  CONTINUE  workspace  is  a  PISETS  workspace  that  has 
been  saved  with  the  >  CONTINUE  command.  This  workspace  is  now  ready  for 
use. 

If  there  is  no  CONTINUE  workspace  in  the  workspace  library  the 
terminal  will  print: 

CLEAR  WS 

This  will  require  that  the  command: 

>  LOAD  PISETS 

be  entered,  following  which  the  terminal  will  print  the  date  and  time 
that  the  PISETS  workspace  was  saved  and  PISETS  will  be  ready  for  use. 

To  quit  either  >  OFF  or  >  CONTINUE  may  be  used. 

If  the  command: 

>  OFF 

is  entered  then  the  present  workspace  contents  will  not  be  saved. 

If  the  command: 

>  CONTINUE 

Is  entered  then  the  present  workspace  contents  will  be  saved  In  the 
CONTINUE  workspace  and  the  previous  contents  of  CONTINUE  will  be  lost. 

0.  USING  THE  WORKSPACE 


The  APL  functions  INPUT,  SET  and  CALC  have  been  defined  In  the 
PISETS  workspace  and  are  used  to  input  the  dimensional  matrix,  choose 


which  sets  of  non -dimensional  parameters  will  be  output  and  to  output 
the  matrix  of  exponents  for  each  of  these  sets. 

The  prompt  >  :  indicates  that  numeric  input  is  required  from  the 
terminal  and  the  prompt  GO:  indicates  that  the  function  has  been 
completed. 

The  function  INPUT  is  used  to  set  up  the  problem.  An  example  using 
the  spring -mass-dashpot  system  developed  earlier  is  shown  on  page  83. 

Entering  the  command: 

INPUT 

activaces  the  function.  In  the  example  INPUT  first  requests  the  number 
of  basic  dimensions  {three)  and  then  the  number  of  dimensioned 
parameters  (seven). 

The  function  then  computes  and  prints  the  maximum  number  of  sets  of 
nondimensi oral  parameters  which  may  exist.  This  is  an  upper  limit  and 
not  all  combinations  will  produce  a  set  of  parameters.  In  the  spring- 
mass-dashpot  example  there  are  35  possible  sets  but  only  24  sets 
actually  exist. 

The  function  now  asks  that  the  dimensional  matrix  be  input.  The 
matrix  is  input  row  by  row,  one  value  at  a  time.  Care  must  be  taken 
when  entering  negative  numbers.  The  negative  sign  is  not  the  same 
symbol  as  the  subtraction  symbol.  The  API  negative  sign  symbol  is  a 
shifted  (uppercase)  2. 

Entering  the  command: 

CALC 

activates  the  APL  function  CALC  which  will  calculate  and  output  an 
exponent  matrix  for  each  set  of  nondimensional  parameters.  Unless  the 
function,  SET,  is  used  innedi ately  prior  to  the  function  CAI.C  all  of  the 
exponents  matrices  that  exist  will  be  output  at  once. 
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The  complete  output  for  all  sets  of  nondimensional  parameters  in 
the  spring-mass-dashpot  example  is  listed  on  pages  86  to  89. 

Entering  command: 

SET 

activates  the  APL  function  SET  which  allows  a  choice  of  which  set  or 
sets  of  nondimensional  parameters  are  output  by  function  CALC. 

The  example  on  page  85  shows  how  SET  is  used  to  choose  one 
particular  set  of  nondimensional  parameters.  The  example  on  page  86 
shows  how  set  ’s  used  to  generate  several  “adjacent"  sets  of  non- 
dimensional  parameters. 

The  words  "PI  TERM  SETS"  refer  to  the  Buckingham  Pi -Theorem  and  are 
synonymous  wit^  "nondimensional  parameter  sets". 

The  user  is  given  the  choice  of  entering  ALL,  1  or  Range. 

If  the  word  ALL  or  the  letter  A  is  entered  the.)  the  next  use  of  CALC 
will  cause  all  the  sets  of  nondimensional  parameters  to  be  output.  This 
Is  identical  to  the  results  of  using  CALC  immediately  after  using  either 
INPUT  or  CALC. 

If  the  number  1  is  entered,  as  shown  in  on  page  84,  then  the  next 
use  of  CALC  will  result  in  the  output  of  only  the  chosen  set  of  non- 
dimensional  parameters.  The  user  chooses  the  nondimensional  parameter 
set  by  choosing  which  of  the  dimensioned  parameters  will  be  considered 
independent.  These  are  the  parameters  which  will  be  fixed  (eliminated) 
when  computing  the  values  In  the  exponent  matrix.  The  user  enters  the 
column  numbers  from  the  dimensional  matrix,  in  ascending  order,  these 
correspond  to  the  chosen  dimensioned  parameters. 

If  the  word  RANGE  or  the  letter  R  is  entered,  as  shown  on  page  85, 
then  the  user  must  choose  a  second  set  of  dimensioned  parameters  to  be 


fixed.  The  next  use  of  CALC  will  result  in  the  output  of  a  range  of 
nondimensi onal  parameter  sets  beginning  with  the  first  set  chosen  and 
ending  with  the  second  set  chosen  and  including  all  sets  "in  between." 
The  numeral  formed  by  the  column  numbers  of  the  second  set  must  be 
"lower"  than  the  numeral  formed  by  the  column  numbers  of  the  first  set. 
for  example  346=  245  >  ?S7  =  2  S  7. 

E .  SAVING  A  WORKSPACE 

For  large  problems  it  may  not  be  desirable  to  print  out  all  possible 
nondimensional  parameter  sets  in  one  sitting.  It  is  also  not  desirable 
to  have  to  input  the  dimensional  matrix  every  time  P I  SETS  is  used. 

When  an  APL  workspace  is  saved  with  the  >  SAVE  conmand  all  of  the 
variables  which  have  been  defined  in  that  workspace  are  saved  along  with 
the  functions.  When  this  workspace  is  later  loaded  with  the  >  LOAO 
command  the  values  of  these  variables  will  be  retained  and  need  not  be 
rei nput. 

After  INPUT  has  been  used  to  set  up  a  problem,  the  workspace  may  be 
saved  with  the  command: 

>  SAVE  PISETS 

and  when  PISETS  is  next  loaded,  work  on  that  problem  can  continue 
without  reuse  of  the  function  INPUT. 

If  it  is  desirable  to  save  the  input  from  more  than  one  problem  then 
the  workspace  for  each  problem  may  be  saved  under  a  different  name.  The 
command : 

>  WSID  NEWNAME 

will  change  the  present  workspace  name  to  NEWNAME.  The  workspace  may 
then  be  saved  with  the  command: 

>  SAVE  NEWNAME 
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When  using  VS  API  under  NERDC  the  >  SAVE  command  will  only  work  if  a 
NEWNAME  workspace  already  exists  in  the  workspace  library.  If  NEWNAME 
is  not  in  the  library  then  the  >  PSAVE  command  must  be  used  in  place  of 
the  >  SAVE  command. 

The  NEWNAME  workspace  may  then  be  substituted  for  the  PISETS 
workspace  when  desired. 

When  using  NERDC,  tne  sign-off  command  >  CONTINUE  causes  the  present 
workspace  to  be  renamed  CONTINUE  and  saved  under  that  name. 

The  command: 

>  LIB 

will  list  the  names  of  all  the  workspaces  which  have  been  saved  in  the 
workspace  library. 

F.  THE  EXPONENT  MATRIX 

The  exponent  matrix,  which  is  obtained  for  each  set  of  non- 
dimensional  parameters,  determines  how  the  dimensioned  parameters  are 
combined  to  form  the  nondimensional  parameters. 

Each  row  of  the  matrix  corresponds  to  a  nondimensional  parameter 
and  each  column  corresponds  to  a  dimensioned  parameter.  Each  non- 
dimensional  parameter  is  formed  by  raising  each  dimensioned  parameter  to 
the  power  shown  in  its  column  and  then  multiplying  all  of  the 
dimensioned  parameters  together. 

Considering  the  spring-mass-dashpot  example  and  the  nondimensional 
parameter  set  that  is  generated  on  page  84.  The  resulting  exponent 
matrix  is  interpreted  as  shown  below: 
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* . I 


y  m  C  k  v0  y0 

PI  1  =  1  0  0  0  0  -1 

PI2  =  0-1  l  o  0  0 

PI3  =  0-1  0  1  0  0 

PI4  8  0  o  0  0  1-1 


0  =  y1VM*k°*vJ*y-1*t°  -  y/y0 

1  *  y(-'*m~1*c 1  Ik^*VQ*y^*t  1  =  ct/m 

2  =  y°*n'1*c1*k°*vj*yj*t2  =  kt2/m 
1  =  y0iin0C0*k0*V(I*y-l*tl  =  v0t/y0 


The  four  terms  PI1,  PI2.  PI3,  and  PI4  are  the  nondimensional  terns 
(pi  terms  j  in  this  set.  The  problem  may  therefore  be  described  by  an 
equation  of  the  form  PI1  =  f(PI2,  P13,  PI4)  or  y/y0  =  f (ct/m,  kt2/m, 
vot/yo) ■ 


i 
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Typical  NERDC  Sign-on  Procedure 


enter  class  nerdc 
class  nerdc  start 

Enter  t  for  TCP,  c  for  CICS,  v  for  vm/cms,  m  for  MUSIC,  1  for  APL 
v 

vm/370  online 

i 


j 


logon  100600843,1 
ENTER  PASSWORD: 
XXXXXXX 


LOGMSG  -  07:10:39  EDT  WEDNESDAY  08/24/83 

*  A  new  release  of  VM/SP  is  now  in  production.  Please  type 

*  "NEW"  for  information  on  changes  in  the  system.  JHB 
LOGAN  AT  16:27:18  EDT  WEDNESDAY  08/24/83 

NERDC  VM/SP3  CMS 


apl 

Z  ( 1A0)  R/0 
X  (1A1)  R/0 
W  (1A2)  R/0 
U  ( 1A3)  R/0 
V  (1A4)  R/0 

R;  T-0. 10/0.21  16:27:26 


APL*  PLUS/1140 
SAVED  16:20:55  08/24/83 


■i 


-  < 


I 

i 

i 

i  i 
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for  T 


Function  INPUT 

Input  for  Spring-Mass-Dashpot  Example 


INPUT 

ENTER  THE  NUMBER  OF  BASIC  DIMENSIONS  USED 
>: 

3 

ENTER  THE  NUMBER  OF  DIMENSIONED  PARAMETERS  USED 
>: 

7 

THERE  ARE  AT  MOST  35  SETS  OF  PI  TERMS 
ENTER  21  VALUES  FOR  THE  DIMENSIONAL  MATRIX  BY  ROUS 
FOR  ROW  1  ENTER  COLUMN  1 
> : 

1 

FOR  ROW  1  ENTER  COLUMN  2 
>: 

0 

FOR  ROW  1  ENTER  COLUMN  3 
>  • 

0 

FOR  ROW  1  ENTER  COLUMN  4 
>: 

0 

FOR  ROW  1  ENTER  COLUMN  5 
> : 

1 

FOR  ROW  1  ENTER  COLUMN  6 
>: 

1 

FOR  ROW  1  ENTER  COLUMN  7 
> : 

0 

FOR  ROW  2  ENTER  COLUMN  1 
> : 

0 

FOR  ROW  2  ENTER  COLUMN  2 
>: 


ETC.  FOR  EACH  VALUE 


FOR  ROW  3  ENTER  COLUMN  6 
> : 

0 

FOR  ROW  3  ENTER  COLUMN  7 
>: 

0 

GO: 
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Function  SET  -  Function  CALC 
Single  Set  of  Nondimensional  Terms 


SET 


HOW  MANY  PI  TERM  SETS  00  YOU  WANT  TO  SEE? 

ENTER  'ALL',  '1',  OR  'RANGE* 

1 

ENTER  THE  COLUMN  NUMBERS  OF  THE  PARAMETERS  TO  BE  FIXED  (ELIMINATED) 
ENTER  3  VALUES  SEPARATED  BY  SPACES 


UNKNOWNS  ELIMINATED  CORRESPOND  TO  DIMENSIONED  PARAMETERS  2  6  7 
THE  MATRIX  OF  EXPONENTS  IS 


10  0  0 
0-110 
0-101 
0  0  0  0 


0-10 
0  0  1 

0  0  2 

1  -1  1 


GO: 
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Function  SET  -  Function  CALC 
Multiple  Sets  of  Nondimensional  Terms 


SET 

HOW  MANY  PI  TERM  SETS  DO  YOU  WANT  TO  SEE? 

ENTER  'ALL',  M  ',  OR  'RANGE' 

RANGE 

ENTER  THE  COLUMN  NUMBERS  OF  THE  PARAMETERS  TO  BE  FIXED  (ELIMINATED) 
ENTER  3  VALUES  SEPARATED  BY  SPACES 
>  : 

3  4  6 

ENTER  THE  COLUMN  NUMBERS  OF  LAST  SET  OF  PARAMETERS 
ENTER  3  VALUES  SEPARATED  BY  SPACES 
> ; 

2  5  7 
GO: 

CALC 

UNKNOWNS  ELIMINATED  CORRESPOND  TO  DIMENSIONED  PARAMETERS  3  4  6 


THE  MATRIX 

OF 

EXPONENT 

c 

IS 

1 

0 

0 

0 

0 

-1 

0 

0 

1 

-2 

1 

0 

0 

0 

0 

0 

1 

-I 

1 

-I 

0 

0 

0 

-1 

1 

c 

0 

1 

UNKNOWNS  ELIMINATED  CORRESPOND  TO  DIMENSIONED  PARAMETERS  345 
THE  MATRIX  OF  EXPONENTS  IS 


1 

i 

0 

-1 

1 

-1 

0 

0 

0 

1 

-2 

1 

0 

0 

0 

0 

0 

-1 

1 

-1 

1 

0 

0 

0 

-1 

1 

0 

0 

1 

UNKNOWNS  ELIMINATED  CORRESPOND  TO  DIMENSIONED  PARAMETERS  2  6  7 
THE  MATRIX  OF  EXPONENTS  IS 


1 

0 

0 

0 

0 

-1 

0 

0 

-1 

1 

0 

0 

0 

1 

0 

-1 

0 

1 

0 

0 

2 

0 

0 

0 

0 

1 

-1 

1 

UNKNOWNS  ELIMINATED  CORRESPOND  TO  DIMENSIONED  PARAMETERS  2  5  7 
THE  MATRIX  OF  EXPONENTS  IS 


1 

0 

0 

0 

-1 

0 

-1 

0 

-1 

1 

0 

0 

0 

1 

0 

-1 

0 

1 

0 

0 

2 

0 

0 

0 

0 

-1 

1 

-1 

GO: 
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Complete  Output  List:  Spring-Mass-Dashpot  Example 


PARAMETERS  5  6  7  FORM  MATRIX  WITH  ZERO  DETERMINANT 


UNKNOWNS  ELIMINATED  CORRESPOND  TO 
THE  MATRIX  OF  EXPONENT  IS 

1  OOOO-IO 
0  10-100-2 
0  0  1  -1  0  0  -1 
0  0  0  0  1  -1  1 

UNKNOWNS  ELIMINATED  CORRESPOND  TO 
THE  MATRIX  OF  EXPONENTS  IS 

1  0  0  0  -1  0  -1 

0  10-100-2 
0  0  1  -1  0  0  -1 

0  0  0  0  -1  1  -1 

UNKNOWNS  ELIMINATED  CORRESPOND  TO 
THE  MATRIX  OF  EXPONENTS  IS 

1  0  0  0  0  -1  0 

0  1  0-1  2-20 

001-11-10 
0  0  0  0  1  -1  1 

UNKNOWNS  ELIMINATED  CORRESPOND  TO 
THE  MATRIX  OF  EXPONENTS  IS 

1  0  0  0  0  -1  0 

0  1  -1  0  0  0  -1 

0  0  -1  1  0  0  1 

0  0  0  0  1  -1  1 

UNKNOWNS  ELIMINATED  CORRESPOND  TO 
THE  MATRIX  OF  EXPONENTS  IS 

1  0  0  0  -1  0  -1 

0  1  -1  0  0  0  -1 

0  0  -1  1  0  0  1 

0  0  0  0  -1  1  -1 

UNKNOWNS  ELIMINATED  CORRESPOND  TO 
THE  MATRIX  OF  EXPONENTS  IS 

1  0  0  0  0  -1  0 

01-101-10 
00-11-110 


DIMENSIONED  PARAMETERS  4  6  7 


DIMENSIONED  PARAMETERS  4  5  7 


DIMENSIONED  PARAMETERS  456 


DIMENSIONED  PARAMETERS  3  6  7 


DIMENSIONED  PARAMETERS  3  5  7 


DIMENSIONED  PARAMETERS  356 


0  0  0 


1  -1  1 


j-j~ 


PARAMETERS  3  4  7  FORM  MATRIX  WITH  ZERO  DETERMINANT 

UNKNOWNS  ELIMINATED  CORRESPOND  TO  DIMENSIONED  PARAMETERS  346 
THE  MATRIX  OF  EXPONENTS  IS 

1  0  0  0  0  -1  0 

0  1-21000 
0  0  1-11-10 
0  0  -1  1  0  0  1 


Complete  Output  List:  Spring-Mass-Dashpot  Example 


UNKNOWNS  ELIMINATED  CORRESPOND  TO  DIMENSIONED  PARAMETERS  3  4  5 
THE  MATRIX  OF  EXPONENTS  IS 

10-11-100 
0  1-21000 
0  0-1  1-1  1  0 
0  0  -1  1  0  0  1 

UNKNOWNS  ELIMINATED  CORRESPOND  TO  DIMENSIONED  PARAMETERS  2  6  7 
THE  MATRIX  OF  EXPONENTS  IS 

1  0  0  0  0  -1  0 

0  -1  1  0  0  0  1 

0-101002 
0  0  0  0  1  -1  1 

UNKNOWNS  ELIMINATED  CORRESPOND  TO  DIMENSIONED  PARAMETERS  2  5  7 
THE  MATRIX  OF  EXPONENTS  IS 

1  0  0  0  -1  0  -1 

0  -1  1  *>  0  0  1 

0-101002 
0000-1  1-1 

UNKNOWNS  ELIMINATED  CORRESPOND  TO  DIMENSIONED  PARAMETERS  256 
THE  MATRIX  OF  EXPONENTS  IS 

1  0  0  0  0  -1  0 

0-1  1  0-1  1  0 

0-101-220 
0  0  0  0  1  -1  1 

PARAMETERS  2  A  7  FORM  MATRIX  WITH  ZERO  DETERMINANT 

UNKNOWNS  ELIMINATED  CORRESPOND  TO  DIMENSIONED  PARAMETERS  246 
THE  MATRIX  OF  EXPONENTS  IS 

2  0  0  0  0  -2  0 

0-12-1000 
0  1  0-1  2-2  0 

0-101002 

'  UNKNOWNS  ELIMINATED  CORRESPOND  TO  DIMENSIONED  PARAMETERS  2  4  5 

THE  MATRIX  OF  EXPONENTS  IS 

2-101-200 

•<  0-12-1000 

0-101-220 

0-101002 

PARAMETERS  2  3  7  FORM  MATRIX  WITH  ZERO  DETERMINANT 

i  UNKNOWNS  ELIMINATED  CORRESPOND  TO  DIMENSIONED  PARAMETERS  236 

THE  MATRIX  OF  EXPONENTS  IS 

!  1  0  0  0  0  -1  0 

0  1-21000 

[  ■  0  1-10  1-10 

I  1  0  -1  1  0  0  0  1 

L 


Complete  Output  List:  Spring-Mass-Dashpot  Example 


UNKNOWNS  ELIMINATED  CORRESPOND  TO  DIMENSIONED  PARAMETERS  235 
THE  MATRIX  OF  EXPONENTS  IS 

1-1  1  0-1  00 

0  1-21000 
0-110-110 
0  -1  1  0  0  0  1 

PARAMETERS  234  FORM  MATRIX  WITH  ZERO  DETERMINANT 
PARAMETERS  1  6  7  FORM  MATRIX  WITH  ZERO  DETERMINANT 


PARAMETERS  1  5  7  FORM  MATRIX  WITH  ZERO  DETERMINANT 

PARAMETERS  1  5  6  FORV  MATRIX  WITH  ZERO  DETERMINANT 

UNKNOWNS  ELIMINATED  CORRESPOND  TO  DIMENSIONED  PARAMETERS  1  4  7 

THE  MATRIX  OF  EXPONENTS  IS 

0  1  0-1  0  0-2 

0  0  1  -1  0  0  -1 

-1  0  0  0  1  0  1 

-1  0  0  0  0  1  0 

PARAMETERS  1  4  6  FORM  MATRIX  WITH  ZERO  DETERMINANT 

UNKNOWNS  ELIMINATED  CORRESPOND  TO  DIMENSIONED  PARAMETERS  1  4  5 

THE  MATRIX  OF  EXPONENTS  IS 

-2  1  0-1  2  0  0 

-101-1100 
-1  0  0  0  0  1  0 

-1  0  0  0  1  0  1 

UNKNOWNS  ELIMINATED  CORRESPOND  TO  DIMENSIONED  PARAMETERS  1  3  7 

THE  MATRIX  OF  EXPONENTS  IS 

0  1  -1  0  0  0  -1 

0  0  -1  1  0  0  1 

-1  0  0  0  1  U  1 

-1  0  0  0  0  1  0 

PARAMETERS  1  3  6  FORM  MATRIX  WITH  ZERO  DETERMINANT 

UNKNOWNS  ELIMINATED  CORRESPOND  TO  DIMENSIONED  PARAMETERS  1  3  5 

THE  MATRIX  OF  EXPONENTS  IS 

-11-10  10  0 

10-11-100 
-1  0  0  0  0  1  0 

-1  0  0  0  1  0  1 

UNKNOWNS  ELIMINATED  CORRESPOND  TO  DIMENSIONED  PARAMETERS  1  3  4 

THE  MATRIX  OF  EXPONENT’S  IS 

0  1-21000 
-10  1-110  0 

-1  0  0  0  0  1  0 

0  0  -1  1  0  0  1 
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Complete  Output  List:  Spri ng-Mass-Dashpot  Example 


UNKNOWNS  ELIMINATED  CORRESPOND  TO  DIMENSIONED  PARAMETERS  1  2  7 

THE  MATRIX  OF  EXPONENTS  IS 

0  -1  1  0  0  0  1 

0-101002 
-1  0  0  0  1  0  1 

-1  0  C  0  0  1  0 

PARAMETERS  1  2  6  FORM  MATRIX  WITH  ZERO  DETERMINANT 

UNKNOWNS  ELIMINATED  CORRESPOND  TO  DIMENSIONED  PARAMETERS  1  2  5 

THE  MATRIX  OF  EXPONENTS  IS 

1-110-100 
2-101-200 
-1  0  0  0  0  1  0 

-1  0  0  0  1  C  1 

UNKNOWNS  ELIMINATED  CORRESPOND  TO  DIMENSIONED  PARAMETERS  1  2  4 

THE  MATRIX  OF  EXPONENTS  IS 

0-12-1000 
-210-1200 
-2  0  0  0  0  2  0 

0-101002 

UNKNOWNS  ELIMINATED  CORRESPOND  TO  DIMENSIONED  PARAMETERS  1  2  3 

THE  MATRIX  OF  EXPONENTS  IS 

0  1-21000 
-1  1-10100 
-1  0  0  0  0  1  0 

0  -1  1  0  0  0  1 

GO: 


PISETS  PROGRAM  LISTING 


Cl  3 

C2D 

133 

C43 

C  5  3 

C63 

C73 

C83 

C  93 

C 103 


vxnput  C133« 

THE  HUMBER  OF  »A«C  DtMEIlSXOMS  USED* 

?eSTSR  THE  HUMBER  OF  DW««0««  PARAMETERS  USED* 

at-Q 

A«.(»»a>f  0 
xa«-»o 

•  THslftE  ARE  ;T"0!Tv;iJ;i'rU  DXMEHSXONAL  MATRZX  ~  ROMS* 

•ENTER  VALUES 


C93  • ENTER  '  » <  T®' 

C 103  ENTSt-0 
C113  XMCROMJ 
C123  ENTD«-EHTr>  +  l 
C133  <(*HT»P)/60 

C 143  entobO 

C153  ihccol; 

C1&3  khtqbemtq+1 


I;::*::-  —  ~ 

C183  ACEHTBf EWTOJpD 
C193  -»<EHTQtQ)/IMCCOL- 

C203  4ZHCROW 

C213  OOJ  *00i 1 

ro23  STOP*-(2*°)~i 

C^33  »asebpowers«-0»4>2*< \°)“1 


C13 
C23 

[33  *  EMTER  •  *  »  1  * 

r  a  *i  CHOICE^  If  Q  __  M4*  •»  b*tdet 

'*•)/«-<-,  ’’I*” 

C&3  -»*NTEn 

\l]  "HVi*  t..  . . .  », »«  ;»«*«-«  TQ  “ 

LOJ  ,  values  seperated  »”  krB 

C93  » ENTER  ’»<?»)»'  VMLUfes 

rifil  powerspD  ___ 

till  »TOP«-+/f'oweHS^0»2*u“row,ERS 

rt*?l  ^{CHOICEb  ’  1  *  ),/*HD  IA»T  SET  or  PARAMETERS' 

S  :bs 

riSl  stoprD 

tl63  tvoPM/2»«-*TOP 

C  17  3  -»EHO 
C 1 B3  R*-1-  * 

ri93  pomeRs«-0  t  *2  *  <  X D)  *  1 

C203  stopp<2*o)-1 
C21  3  eHD  • '  GO  *  * 


VSET  C  0  3  v 

V  SET  TERM  SETS  DO  TOU  WANT  TO  SEE?' 

•HOW  MAIU  PZ  TERM  St 

ENTER  :  ...  ......  OR  *  '  RANGE’  ’  * 


<  ALL  • 
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LISTING  CONTINUED 


VDZT  £Q3v 
y  Zf-DET  F  ;  v  ;  N  ;  J 

[13 

[23  v«-o »  ( >■’- 1 )  f  1 

£33  I  FC  » 1 3 )  l  ["  /  I  FC  » 1  3 

143  -»<l=->)/7 

£53  PCI ,  j;  ]FF  ,:j,  1  ;3 

C63  =  «-- 2 

C73  4(iE-io< |PC1»13)/10 
C83  =*-0 

£93  ^RETURN 

CJ03  2«.2xFClfl3 
[113  PCl?3«-PCi;3^PCi;i3 
C  1  2  ]  p<-p-Pl!13*.xFC1»3 

C13  3  PW/C13  v/f 
C  1 4  3  -»(0<,»v*--14v)/3 

CIS  3  xeturhj 
v 


1 


I 


I 

I 


yC«L.C  COUv 
y  CALC 

£13  ehcope;  1  • 

[23  ^((+/F-OWtF:r.)>STOP)/EMD 

£33  I»K  (o,'2)t  +  /foweP5)/xq 

C 4 3  -*<  (  |Zrdf.t  ACfii' 3)<0.0001)/2ERo 

£f,3  1  UMKHO'JMS  cLIMIHarES  CORRESFONH  TO  DIMENSIONED  PARAMETERS  ’  r(TIJ) 

£63  HiSK-iOfiD  )/iq 
T73  Ef(Q,RI)rO 

C  8  3  ECIE>;  3  m  <  3qC  ;  103)  +  .  x(-A£;wi&3) 

£93  ECNir.;3<-<FXfPi)fip»>lfO 

£103  ’THE  MATRIX  GF  EXPONENTS  IS  1 

[11]  $ ( ( I E) >0 .01 >  x£x  I  2 

C 123  ■♦next 

£133  zeroj 

£143  'PARAMETERS  ',(TXD),'  r'J-M  MATRIX  WITH  ZERO  DETERMINANT* 

£153  next; 

£163  ned+1 

£173  LO '■p; 

C183  POWERS  £N34-2XPOWERS£N3 

C193  ■*  <  POWERS  £H]^POWERS£N-  13  )  /  ENCODE 

£203  POWEXSfH](.»ASE[H3 

£213  hj-n-1 

£223  -+UO  OR 

£233  end; -  go;  • 

C243  POWERS«.PAS£ 

£25  3  stofm2,c,)-i 

V 


- «. 
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